Precise measurement of the top quark mass in the lepton+jets topology at CDF II 
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955 pb -1 of integrated luminosity, we achieve the single most precise measurement of the top quark 
mass, 170.8 ± 2.2 (stat.) ± 1.4 (syst.) GeV/c 2 . 
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The top quark is the heaviest known elementary parti- 
cle, about 40 times more massive than the next-heaviest 
fcrmion. Its large mass plays an important role in loop 
corrections to several electroweak observables predicted 
by the standard model. In conjunction with measure- 
ments of the W boson mass, a precise measurement of 
the top quark mass, m t , constrains the mass of the as 
yet unobserved Higgs boson []]]. Moreover, by compar- 
ing precision electroweak measurements to predictions in- 
cluding the relevant loop corrections, a precise measure- 
ment of m t can help constrain contributions from physics 
beyond the standard model. Measuring m t to the high- 
est achievable precision is therefore one of the main goals 
of the experiments operating at the Fermilab Tevatron 
collider. 

At the Tevatron, top quarks are mainly produced in 
pairs. They decay essentially 100% of the time into a W 
boson and a b quark, with the W decaying into quarks or 
leptons. The result presented here uses the lepton+jets 
channel, where one W decays into two quarks, and the 
other decays to an electron or a muon and the corre- 
sponding neutrino. In the past, this channel has pro- 
vided the most precise measurements of mt, and recent 
measurements can be found in Ref. Q. 

In this letter we report the single most precise mea- 
surement of the top quark mass from Tevatron proton- 
antiproton collisions at -^3=1.96 TeV, using 955 pb -1 of 
integrated luminosity collected with the CDF II detector 
from February 2002 to March 2006. The CDF II detector 
is a general-purpose particle detector and is described in 
detail elsewhere [3]. It has a solenoidal charged particle 
spectrometer, consisting of 7-8 layers of silicon microstrip 
detectors and a cylindrical drift chamber immersed in 
a 1.4 T magnetic field, a segmented sampling calorime- 
ter, and a set of charged particle detectors outside the 
calorimeter used to identify muon candidates. We use a 
right-handed cylindrical coordinate system with the ori- 
gin in the center of the detector, where 9 and (j> are the 
polar and azimuthal angles, respectively, and pseudora- 
pidity is defined as rj = — In tan(#/2). Transverse energy 
and momentum are Et=E sin 6 and px—psmO, respec- 
tively, where E and p are energy and momentum. 

Events in the lepton+jets decay channel are selected 
to have a single, isolated electron or muon candidate 
with large transverse energy, large imbalance in trans- 
verse momentum in the event (missing transverse en- 
ergy i$x) as expected from the undetectable neu- 
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trino, and exactly four jets with large transverse energy. 
Jets are reconstructed using a cone algorithm with radius 
AR = ^(Ar)) 2 + (A<f)) 2 = 0.4. Of these jets, we require 
at least one to have originated from a b quark by using an 
algorithm that identifies a long-lived B hadron through 
the presence of a displaced vertex (b tag) Back- 
grounds to the it signal arise from multi-jet QCD pro- 
duction (non-W), W production in association with jets 
(W+jets), and electroweak backgrounds (EWK) com- 
posed of diboson (WW, WZ, ZZ) and single top pro- 
duction. VF+jets background events include jets with 
real b flavor as well as light flavor jets incorrectly identi- 
fied as b jets. To remove the non-14 7 backgrounds, where 
$t is due to mis-measured jet energies, we require j^r not 
to be aligned with the highest energy jet by a suitable 
requirement on A(f> between this jet and $r- Table U 
summarizes the selection criteria used in this analysis, 
and a more detailed description can be found in Ref. [|| . 
We select 167 events of which we expect about 85% to 
be it events. Table Ql] shows the expected sample com- 
position determined with 318 pb -1 [6(, scaled to 955 pb -1 
and assuming a it cross section of 8.0 pb. Detailed back- 
ground estimate studies justify this scaling, and residual 
differences are absorbed in the systematic background 
uncertainty. 

TABLE I: Event selection criteria 0. 



Lepton Et >20 GeV (electron, muon), \rj\ < 1 

Jets exactly 4 with E T > 15 GeV, \r)\ < 2.0 
# T > 20 GeV, calculated over \r)\ < 3.6 

6 tag jets > 1 from a secondary vertex, \r]\ < 1.5 



Non-VK veto 0.5 < A<j> < 2.5 when # T < 30 GeV 



TABLE II: Background composition and expected number of 
tt candidates. All uncertainties are statistical only. 



Source 


Expected number of events 


W+jets 


14.5 ± 5.1 


non-VK 


5.2 ± 2.6 


EWK 


2.2 ± 0.5 


Total 


22.0 ± 8.2 


it (cr=8.0 pb, m t =170 GeV/c 2 ) 


145.1 ± 16.5 


Data 


167 



We analyze the selected events using a likelihood tech- 
nique that relies on calculations of probability densities 
based on matrix elements for the signal (tt) and dominant 
background (W+jets) processes [7J. The backgrounds 
other than VF+jets are found to be adequately described 
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by the W+jets probability density. Given a set of ob- 
served variables, x, and underlying partonic quantities, 
y, the signal and background probability densities are 
constructed by integrating over the appropriate parton- 
level differential cross section, da(y)/dy, convolved with 
parton distribution functions (PDFs) and detector reso- 
lution effects: 



P(x) 



E 



jet perm. 



do{y) 
dy 



f(qi)f(Q2)dqidq 2 W(x, y)dy. 



(1) 

The PDFs (f{qi) and /(tft)) take into account the fla- 
vors of colliding quark and antiquark and are given by 
CTEQ5L [8]. The detector resolution effects are de- 
scribed by a transfer function W{x 1 y) relating x to y. 
The momenta of the leptons and the angles of jets and 
leptons are taken to be exactly measured, and therefore 
W(x,y) for these quantities is given by the product of 
Dirac delta functions. The non-trivial part of W(x,y) 
maps parton energies to measured jet energies after cor- 
rection for instrumental detector effects [9(. This map- 
ping is obtained by parameterizing the jet response in 
fully simulated ti events created by the Monte Carlo 
(MC) generator pythia [13], and including the effects 
of radiation, hadronization, measurement resolution, and 
energy omitted from the jet cone by the reconstruction 
algorithm. The ti and VF+jets probability densities, P t i 
and Pw+jets , include all possible permutations of match- 
ing jets with partons as well as all possible longitudinal 
momenta for the neutrino in the W decay. The permu- 
tations are reduced to six or two by exploiting b tagging 
information (single-tag or double-tag, respectively). We 
use different transfer functions for light quark jets and b 
jets, depending on the flavor of the parton assigned to the 
jet. In calculating da(y)/dy, P t i uses the leading order 
matrix element of the qq — > ti process and Pw+jets 
uses the sum of matrix elements of the W+A jets sub- 
routines encoded in the vecbos Monte Carlo generator 

The final state described by da(y)/dy contains 6 par- 
ticles, which introduces 20 integration variables in Eq. [TJ 
including the longitudinal momenta of the incoming 
quarks. By imposing energy-momentum conservation, in 
conjunction with the Dirac delta functions in W{x 1 y) we 
reduce the dimensionality of the remaining integration to 
five. The integration in Pw+jets is performed over the 
energies of the outgoing partons and the invariant mass of 
the leptonically decaying W using a Monte Carlo tech- 
nique. In order to reduce the calculation time for P t t, 
we integrate over the following variables: the invariant 
masses of i, i, W + , and W~ , and the energy of one of 
the quarks from the hadronic W decay. Our method 
includes two additional integrations over the transverse 
momentum components of the ti system. The integration 
in P t i uses the numerical integration code VEGAS 131 ]. 

The largest potential systematic uncertainty in this 



measurement arises from the energy scale of jets. To 
decrease this uncertainty, we exploit the fact that the 
hadronically decaying W provides an in situ constraint 
of the jet energy scale, as the two jets should form an 
invariant mass consistent with the precisely known mass 
of the W boson Q. The jet energy scale and the mass of 
the top quark are simultaneously determined from a two- 
dimensional likelihood that includes their correlation. A 
salient feature of this method is that the uncertainty due 
to the jet energy scale will be reduced with increasing 
statistics. Thus P t i is evaluated as a function of mt and 
an assumed jet energy scale factor Jjes = Ejet/Ejet, 
where is the observed jet energy and Ej et is the 

true jet energy. 

To extract m t and fjES from the data, we build a 
likelihood function for N selected events by adding P t t 
and Pw+jets for each event. The combined likelihood 
is minimized with respect to three variables: mt, fjES, 
and C s , the fraction of events consistent with our ti signal 
hypothesis. The likelihood for TV events is given by 

£(xi,x 2 , ...,x N ;m t , fjES,C s ) = 



-N(C s (A t{ (m t ,fj ES )) + (l-C s )(A w+jBts (fjEs))) 



N 



JJ[C s P t i(x; m t , fjEs) + (1 - C s )P w+iets (x; f JE s) 



(2) 



where the first factor arises from the Poisson extension of 
the likelihood and normalizes the combined event prob- 
ability density, and (^4) refers to the mean acceptance 
for ti or W^+jets events. We use fully simulated MC ti 
and ly+jets events to determine the functional form of 
(A). Pw+jets is evaluated at the central jet energy scale 
factor, /,/£s=l- The f JE s dependence of P w + 3 ets is de- 
termined by varying the input Jjes m MC event samples 
(fjEs) ano - ^ parameterizing the average likelihood re- 
sponse as a function of fjES- We use the rrtt dependence 
of the theoretical leading order ti cross section to normal- 
ize P t i- Because we use a leading order matrix element 
to calculate P t i , we find that ti events where at least one 
of the four reconstructed jets cannot be matched to a 
parton from the ti decay within AR < 0.4 behave like 
background events. As a consequence, a pure sample of 
ti events yields C s of 0.8. The quoted C s values are cor- 
rected for this effect. For each event P t i is evaluated in 
increments of 2 GeV/c 2 in m t and 0.02 in Jjes- At each 
point of this grid we fit the entire sample of N events 
according to Eq. [2] and the most likely value of C s is 
determined using MINUIT [3]. The optimal parameters 
m t and Jjes are obtained by fitting the likelihood using 
a two-dimensional Gaussian. The statistical uncertainty 
on mt includes the uncertainty on fjES- 

The performance of the analysis is tested by extract- 
ing mt from MC pseudo-experiments containing ti sig- 
nal samples with various input top quark masses (m^ ) 
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and background samples described in Table II. The signal 
and electroweak background samples are generated using 
herwig [HI. The IF+jets background is generated using 
ALPGEN [101 ] with hadronization and fragmentation done 
by herwig. The non-VT" background is extracted from 
an independent data sample. All of the MC samples are 
processed by the CDF detector simulation. We construct 
pseudo-experiments of signal and background events by 
fluctuating the number of events around the values shown 
in Tabic [ill Figure QJa) shows that the fitted Gaus- 
sian mean m t extracted from 200 pseudo-experiments 
per point is unbiased with respect to mf' IC up to the 
statistical uncertainty of 0.21 GeV/c 2 , which is taken as 
a systematic uncertainty. Similar tests are performed for 
the output of fjES- In this case, we find that a bias of 
+4% in Jjes is present, independent of m t . We correct 
for this bias to properly interpret the output of fjES- 
Figure QJb) shows the top mass pull width, defined as 
the Gaussian a of the top mass residual (m t — mf' IC ) di- 
vided by the uncertainty in each pseudo-experiment a„ H , 
as a function of m\ lc . The pull width is 3%±2% larger 
than one on average, and thus, the statistical uncertainty 
is scaled up by 3%. 



ti events and is consistent with the expectation shown in 
Table HI1 Monte Carlo tests have shown that the resulting 
mt is stable over a wide range of sample purities. Fig- 
ure [3] shows comparisons of two representative kinematic 
quantities between data and simulation using /j|j^=0.99 
and mf c =170 GeV/c 2 . 
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FIG. 2: Contours of likelihood evaluated over the 955 pb" 
event sample. Our measurement is indicated by the X. 




165 170 175 180 165 170 175 180 

m t MC (GeV/c 2 ) m, MC (GeV/c 2 ) 

FIG. 1: Results of pseudo-experiment tests, (a) Difference 
between the measured m t and the input top quark mass in the 
MC event sample (mf IC ), as a function of mt' c . (b) Gaussian 
a of pull distributions (see text), as a function of m^ 10 . The 
plots include results using MC event samples with different 
values of fjES (fjEs)- The weighted average po is indicated 
by the solid horizontal line. The dashed line indicates an 
example of an unbiased result. 

Applying this method to data, we measure the top 
quark mass to be m t =170.8±2.2(stat.), and the Jjes 
scale to be /j_ES=0.99±0.02(stat.) in good agreement 
with the reference scale from the default CDF calibra- 
tion Figure [2] shows the fitted two-dimensional likeli- 
hood with AlnL contours. The statistical uncertainty is 
taken from the maximum and minimum m t values on the 
A In £=0.5 contour. We find a correlation coefficient of 
0.32 between mt and /jes- The fit yields a signal frac- 
tion C s =0.84±0.10(stat.), which corresponds to 140±17 




m t (GeV/c 2 ) two-jet mass (GeV/c 2 ) 



FIG. 3: Comparison of two kinematic variables for data and 
MC using ff E c s =0.99 and mf c =170 GeV/c 2 . (a) Most prob- 
able value of mt for each event extracted from evaluating P t j 
at fjES~l- (b) Invariant mass of the pair of jets assigned 
as W decay products calculated using the most probable per- 
mutation at the most probable value of mt and fjES in each 
event evaluated from P t j. The backgrounds contain all con- 
tributions shown in Table ITT1 

The sources of systematic uncertainty are listed in Ta- 
ble IIIII To first order, Jjes is already included in the 
statistical uncertainty, but we also consider a dependence 
of Jjes on the px and 77 of the jets (residual jet energy 
scale) using the dependence found in other studies @]. 
Another uncertainty is calculated from possible differ- 
ences in Jjes between b jets and light quark jets. The 
generator uncertainty takes into account differences in 
parton showering and jet fragmentation between two dif- 
ferent MC programs used to generate ti events, pythia 
and herwig. Variations in initial- and final-state radia- 
tion (ISR, FSR), constrained by studies using Drell-Yan 
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data, are also considered [2J . PDF uncertainties are eval- 
uated using MC samples generated with MRST [l5[ and 
the full set of eigenvectors from CTEQ6M Q. Systematic 
effects on the modeling of the background samples in- 
clude fluctuations in the total background contribution, 
relative contributions from individual background pro- 
cesses, and variations due to the Q 2 scale used in W+jets 
simulation. We also include effects from the uncertainties 
in the simulated lepton px, the dependence of b tagging 
with jet pr, and the effect of the limited MC event sam- 
ples used in the analysis. Finally, we include possible 
mis-modeling of multiple interactions in the simulation 
at high luminosity. The sum in quadrature of all system- 
atic uncertainties is 1.4 GeV/c 2 . 
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TABLE III: Summary of systematic uncertainties. 



Source 


Uncertainty (GeV/c 2 ) 


Residual jet energy scale 


0.4 


6 jet energy scale 


0.6 


Generator 


0.2 


ISR/FSR 


1.1 


PDFs 


0.1 


Background 


0.2 


Lepton pr 


0.2 


b tag pr dependence 


0.3 


Monte Carlo statistics 


0.2 


Multiple interactions 


0.1 


Total 


1.4 



In summary, we present a measurement of the top 
quark mass in the lepton+jets channel using 955 pb" 1 
of data collected by the CDF experiment. A matrix ele- 
ment analysis was used with an in situ measurement of 
the jet energy scale. We measure 

rat = 170.8 ± 2.2 (stat.) ± 1.4 (syst.) GeV/c 2 (3) 

where the statistical uncertainty includes the uncertainty 
of 1.5 GeV/c 2 due to the jet energy scale. With a total 
uncertainty of 1.5%, this result is the most precise mea- 
surement of the top quark mass to date and is a 35% 
improvement over the previous best measurement 0]. 
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